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DYNAMICAL SYSTEMS GENERATED BY A GONOSOMAL 
EVOLUTION OPERATOR 

U. A. ROZIKOV, R. VARRO 


Abstract. In this paper we consider discrete-time dynamical systems generated by 
gonosomal evolution operators of sex linked inheritance. Mainly we study dynamical 
systems of a hemophilia, which biologically is a group of hereditary genetic disorders 
that impair the body’s ability to control blood clotting or coagulation, which is 
used to stop bleeding when a blood vessel is broken. We give an algebraic model 
of the biological system corresponding to the hemophilia. The evolution of such 
system is studied by a nonlinear (quadratic) gonosomal operator. In a general setting, 
this operator is considered as a mapping from R", n > 2 to itself. In particular, 
for a gonosomal operator at n = 4 we explicitly give all (two) fixed points. Then 
limit points of the trajectories of the corresponding dynamical system are studied. 
Moreover we consider a normalized version of the gonosomal operator. In the case 
n = 4, for the normalized gonosomal operator we show uniqueness of fixed point and 
study limit points of the dynamical system. 
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1. Introduction. 

In biology it is important a proper understanding of living populations at all levels. 
The relevant mathematics undoubtedly requires of nonlinear analysis, in particular 
a nonlinear dynamical system, compounded stochastic processes modeling, and the 
creative implementation of the computer methodology. 

The action of genes is manifested statistically in sufficiently large communities of 
matching individuals (belonging to the same species). These communities are called 
populations [8]. The population exists not only in space but also in time, i.e. it has its 
own life cycle. The basis for this phenomenon is reproduction by mating. Mating in a 
population can be free or subject to certain restrictions. 

The whole population in space and time comprises discrete generations Fq,F\,... 
The generation F n+ ± is the set of individuals whose parents belong to the F n genera¬ 
tion. A state of a population is a distribution of probabilities of the different types of 
organisms in every generation. 

A type partition is called differentiation. The simplest example is sex differentiation. 
In bisexual population any kind of differentiation must agree with the sex differentia¬ 
tion, i.e. all the organisms of one type must belong to the same sex. Thus, it is possible 
to speak of male and female types (see for example EL EL P20 for mathematical models 
of bisexual population). 
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In many biological systems, sex is determined genetically: males and females have 
different alleles or even different genes that specify their sexual morphology. In animals, 
this is often accompanied by chromosomal differences. Determination genetically is 
generally through chromosome combinations of XY (for example: humans, mammals), 
ZW (birds), XO (in this variant of the XY system, females have two copies of the 
sex chromosome (XX) but males have only one (ATO). The 0 denotes the absence 
of a second sex chromosome. There are some sex linked systems which depends on 
temperature and even some of systems have sex change phenomenon (see [TO] for a 
detailed review.) We note that the behavior of sex-linked system can be investigated 
by studying of nonlinear dynamical systems, such systems are not fully understood yet. 
A search in MathSciNet gives about 15 mathematical papers which are related to sex 
linked models (see for example, m, 0 , 0 , m)- In papers m , m we attempted to 
introduce thermodynamic methods in biology. In [TO] an algebra associated to a sex 
change is constructed. 

In this paper we consider evolution (dynamical system) of a hemophilia. Recall that 
hemophilia is a group of hereditary genetic disorders that impair the body’s ability 
to control blood clotting or coagulation, which is used to stop bleeding when a blood 
vessel is broken. 

In the next section we give a mathematical model of the biological system corre¬ 
sponding to the hemophilia. The evolution of such system will be given by a nonlinear 
(quadratic) evolution operator which is called a gonosomal operator. Thus study of 
the biological system is reduced to the study of the nonlinear dynamical system gen¬ 
erated by the gonosomal operator. In a general setting, this operator is considered as 
a mapping from M n , n > 2 to itself. In Section 4 we give some detailed properties of 
the dynamical system. In particular, for a gonosomal operator at n = 4 we explicitly 
give all (two) fixed points. Then limit points of the trajectories of the corresponding 
dynamical system are studied. In the last section we consider the normalized version 
of the gonosomal operator. In the case n = 4, for the normalized gonosomal operator 
we show uniqueness of fixed point and study limit points of the dynamical system. 


2. Bisexual population: Gonosomal evolution operator 

Type partition is called differentiation. The simplest example is sex differentiation. 
In bisexual population (BP) any kind of differentiation must agree with the sex differ¬ 
entiation, i.e. all the organisms of one type must belong to the same sex. Thus, it is 
possible to speak of male and female types. 

In many cases, the sex determination is genetic, in particular, it is controlled by two 
chromosomes called gonosomes. Gonosomal inheritance is a mode of inheritance that 
is observed for traits related to a gene encoded on the sex chromosomes. 

Let us discuss one example of sex-linked inheritance. Haemophilia is a lethal recessive 
A-linked disorder: a female carrying two alleles for hemophilia die. Therefore if we 
denote by X h the gonosome X carrying the hemophilia, there are only two female 
genotypes: XX and XX h ( X h X h is lethal) and two male genotypes: XY and X h Y. 
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We have four types of crosses: 

XX x XY >-> iXX, iXX; 

XX x X h Y >-» ±XX h , \XY- 

2 2 ( 2 . 1 ) 

xx h xxfm ±xx, \xx h , ±xy, ±x h y; 

xx h x x h y >-► ixx' 1 , ±xy, ±x h y. 

Let F = {XX,XX fi j and M = {XL, X h y} be sets of genotypes. Assume state of the 
set F is given by a real vector (x,y) and state of M by a real vector (u,v). Then a 
state of FUM is given by the vector s = ( x , y, u, v ) £ R 4 . If s' = (x', y', u', v') is a state 
of the system F U M in the next generation then by the rule (12.II) we get the evolution 
operator W : R 4 —> R 4 defined by 


' x' = \xu + \yu 
y' = ±xv + \yu + \yv 
u' = \xu + \xv + \yu + |yv 
. v' = \yu + \yv. 


( 2 . 2 ) 


This example can be generalized: suppose that the set of female types is F = 
{1,2,..., n} and the set of male types is M = {1, 2,..., is}. Let x = (xi ,..., x n ) £ R n 
be a state of F and y = (yi,..., y u ) £ IT be a state of M. 

Consider 7 and 7 ^ as some inheritance real coefficients (not necessary probabil¬ 
ities) with 

n v 

E4S+E^5 ) = 1 - (2-3) 

3 = 1 i=i 

Consider an evolution operator W : R n+iy —> R n+l ' defined as 


W : 


Xj = } 2 i,k= 1 Tik 'jXiyk, J = 1, • • •, n 

f \~~yTi.is (m) 7 1 

Vl = Di,fc=1 TikJXiVk, 1 = 1, • • • , P. 


(2.4) 


This operator is called gonosomal evolution operator. This means that the association 
s = (x,y) £ R' 1 " 1 "" —> s' = (x',y r ) £ R n+l/ defines a map W. The population evolves 
by starting from an arbitrary state s, then passing to the state s' = W(s) (in the 
next ’generation’), then to the state s" = W(W(s)), and so on. Thus, states of the 
population described by the following discrete-time dynamical system 


s (0) , s (1) = IT(s (0) ), s (2) = IT 2 (s (0) ), s (3) = W 3 (s (0) ), 


(2.5) 


where s® £ R n+ " is a given initial point and W n (s) = W(W(...W(s))...) denotes the 

n 

n times iteration of W to s. 

The main problem for a given dynamical system is to describe the limit points of 
the trajectory for arbitrary given s®. 
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3. Dynamical system generated by the operator (12.2|) 

Note that operator (|2.4)l describes evolution of a hemophilia. The dynamical systems 
generated by gonosomal operator (12.41) is complicated. In this paper we study the 
dynamical system generated by gonosomal operator (12.21) . which is a particular case of 
(12.41) . obtained by n = u = 2 and the following coefficients: 


7& = \ 
112,1 = 0 
i2i\ = i 

7$1 = 0 


111,2 = 0 

7l 2,2 = 5 
Jf) - I 

721.2 — 4 

-(/) - I 

722.2 — 3 


(m) 

7n,i 

(m) 

712.1 

(m) 

721.1 

(m) 

722.1 


33 = o 
33 = 0 

7(3 = 3 

7(3 = 3 


3.1. Fixed points. 

A point s is called fixed point if W(s) = s. Let us find all fixed points of W given 
by (12.21) . i.e. we solve the following system of equations 


W : 


x = \xu + \yu 
y = \xv + \yu + \yv 
u = \xu + \xv + \yu + \yv 


(3.1) 


2 Hi | 2 " 

V = \yu+ \yv. 

First it is easy to see that so = (0,0, 0,0) is a solution to system (13.11) . 

To find another solution from the first equation of this system we get (4 — 2 u)x = yu. 
Assuming u = 2 from this equation we get y = 0 then the last equation of (13.11) gives 
v = 0 consequently from the third equation of the system we get x = 2. Thus we 
obtained the solution S 2 = (2,0,2,0). 

From the last equation we get (12 — 4 y)v = 3 yu, assume first that y = 3 then this 
equation gives u = 0, consequently the first equation of the system (13.ip gives x = 0. 
Then from the second equation we get v = 3. But these values do not satisfy the third 
equation of the system. Assume now u ^ 2 and y / 3 then we have 

yu 3 yu 


x = 


v = 


(3.2) 


4 — 2ri 12-4 y 

Using (13.21) from the second and third equations of the system (13.11) we obtain 

J 16(3 — y)(2 — u) = 3rt(8 — 4 u + yu) 

{ 16(3-y)(2-u) =y(24-yu). 

From the first equation of the last system we find 

_ 12(rt 2 - 6 u + 8) 

V ~ 3 u 2 - 16m + 32 ' 

Substituting this to the second equation of the last system we obtain the following 
equation 

(u - 2) 2 (u - 8) (3 u 2 - 14u + 24) = 0. 
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This equation gives u = 2, u = 8 and 3 u 2 — 14ii + 24 = 0. For case u = 2 we have 
solution S 2 mentioned above. The case u = 8 gives y = 3 which does not give solution 
of the system (13.111 as was discussed above. Thus only remains 3 u 2 — 14 u + 24 = 0 
which does not have real solutions. 

We proved the following 

Proposition 1. The gonosomal operator 12.2 I) has exactly two fixed points: sq = 
(0,0, 0,0) and S 2 = (2,0, 2,0). 

3.2. The type of the fixed points. 

Now we shall examine the type of the fixed points. 

Definition 1. (see m- A fixed point s of the operator W is called hyperbolic if its 
Jacobian J at s has no eigenvalues on the unit circle. 

Definition 2. (see m- A hyperbolic fixed point s is called: 

• attracting if all the eigenvalues of the Jacobi matrix J(s) are less than 1 in 
absolute value; 

• repelling if all the eigenvalues of the Jacobi matrix J(s) are greater than 1 in 
absolute value; 

• a saddle otherwise. 


To find the type of a fixed point of the operator (12.2jl we write the Jacobi matrix: 

( \u \u + 0 \ 


J(s) = Jw = 


V 


\u 

\ u 

\ x + \y 

0 

h v 

{u + ±v 

\ y 

h x + b 

\ u + \ v 

h + b 

\ x + \y 

\ x + b 

0 

\u + \v 

\y 

b 




It is easy to see that J(sq) has all eigenvalues equal to 0, therefore so is an attracting 
point. 

The Jacobian J(s 2 ) has eigenvalues — ^,0,1,2, therefore the fixed point is not hy¬ 
perbolic. 


3.3. Dynamics on invariant sets. 

A set A is called invariant with respect to IT if IT (A) C A. 

Denote 

O = {(0,0, u, v ) € M' 1 : u, v € R} U {(x, y, 0, 0) € M 4 : x, y € M}, 
I = {s = (x, y, u, v) € M 4 : y = v = 0}, 

J = {s € / : x = u}, 

P = {s = (x, y, u, v) € M 4 : x > 0, y > 0, u > 0, v > 0}, 

Q a = {s = (x,y,u,v) € P : x + y + u + v < a}, a£ [0,4], 

A f = {s = (x, y, u, v) € M 4 : x < 0, y < 0, u < 0, v < 0}, 

No = {s = (x, y, u, v) € M 4 : x < 0, y < 0, u > 0, v > 0}, 

Ni = {s = (x, y, u, v) € M 4 : x > 0, y > 0, u < 0, v < 0}. 
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Lemma 1.1) The sets I, J, P and Q a (a G [0,4],) are invariant with respect to W. 

2) W{0) = {( 0,0,0,0)}. 

3) W(Q a ) C Q a 2 /4 . 

4) W(AT) C P. 

5) W{Af 0 ) c Af, W{Af 1 ) c Af. 


Proof. We give the proof for Q a , for other sets it simply follows from ([2.2D . Take any 
s = (x , y, u, v ) G Q a then we have 0 < x + y < a and 0 < u + v < a. From (12.21) we get 
x' > 0, y' > 0, v! > 0, v 1 > 0 and 

( \ 2 9 

x+y+u+v\ a 

-2-) - T' 

Thus s' = ( x ', y ', u', u') G C Q a - □ 


Reduce W on J then we get the mapping x' = f{x) = ^x 2 . This function has two 
fixed points x = 0 and x = 2. Moreover, 0 is attractive (/'(0) = 0 < 1) and 2 is repeller 
(/'(2) = 2 > 1). Take an initial point Xq € J and iterate the function /, then we get 


x n = f n (x 0 ) 


o—(l+2+2 2 ...+2 n_1 ) 2 71 

2 X 0 


2 n +l 2 71 _ 


Xn 


= 2 



Hence we have 


lim x n = < 

n—>oo 


Now reduce the operator W on /: 


V : 


0 , if |x 0 | < 2 
2 , if |x 0 | = 2 

+oo, if |x 0 | > 2 . 

x' = \xu 
v! = \xu. 


Thus for any to = (xo,0, uq,0) G I we have W(xo, 0, no, 0) G J. Consequently, we have 
full characterization of the dynamical system on the invariant set J, i.e., we proved the 
following 


Proposition 2. For any initial point Q = (xo, 0, uo, 0) £ I we have 

' ( 0 , 0 , 0 , 0 ) if |x 0 u 0 | < 4 


lim W n (t 0 ) = 

n—too 


( 2 , 0 , 2 , 0 ) if |x 0 uo| = 4 

Too, if |x 0 uo| > 4. 


Let us now consider the dynamical system on the other sets (which may intersect 
with J). 

Lemma 2. Let a G [0,4). Then for any initial point s = (x,y,u,v) G Q a we. have 


lim W n (s) = (0,0, 0,0). 

n—>oo 


(3.3) 
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Proof. Let /(a) = a 2 /4. By LemmaUwe have 

W n {Q a ) C W n - l (Q f[a) ) C W n ~ 2 (Q f 2 {a) ) C • • • C Q/n^). 

It is easy to see that f(x) has two fixed points 0 and 4. Moreover, 0 is attracting point 
and 4 is repelling point. For any a G [0,4) we have lim n f n (a) = 0. Consequently, 
we get Hindoo W n (Q a ) C Q 0 = {(0,0, 0,0)}. □ 

Lemma 3. For an initial point s = (x, y, u, v ) € Qa the following hold 

i. if there is k > 0 such that fz 0 then \3.3f) is satisfied, 

ii. if y( fc )y( fe ) = 0 for any k > 0 then 

lim W n (s) = (2,0,2,0), (3.4) 

n—>oo 

where y^ and v^ are second and fourth coordinates of the vector W k (s). 

Proof. In the case a = 4 we have 0<x + y + u + u<4. From this inequality it follows 
that 

0 < x + y < 2 or 0 < u + v < 2, 

if both x + y and u + v large than 2 then their sum is large than 4. Without loss of 
generality we assume that t = x + y < 2 then u + v < 4 — t. Hence we have 

f 4, if t = 2 

x' + y' + v! + v' = (x + y)(u + v) < t{4 — t) = < 

{ < 4, if i < 2 

hence W(s) G Qt( 4 -t), i- e -i for t < 2 the case is reduced to the case of a < 4. Consider 
now the case t = 2. Then if u + v < 2 we can reduce the case to the case a < 4. But if 

t = x + y = 2 and u + v = 2 (3.5) 

then from (12.21) we get 

x' + y' + v! + v' = 4, x' + y' = 2 — — and u' + u / = 2 + —. (3.6) 

6 6 

Consequently, if yv / 0 then W 2 (s) G Q 4 _(ee_) 2 . By (13.51) we have 0 < y < 2 and 

0 < v < 2, hence 0 < 4 — (^) 2 < 4. Thus condition (13.51) together with yv / 0, by 

Lemma [2] gives (13.31) . 

Let now yv = 0 then (13.61) is reduced to the case (13.51 ). Repeating above argument 
we see that if y'v' / 0 then VF 3 (s) G otherwise we iterate the argument 

again. By this way one can show that if (13.51) is satisfied and there exists k > 0 such 
that 'i/^viF fz 0 then we have (13.31) . 

Suppose now (13.51) is satisfied and 

y{k) v (k) _ Q £ or an y fc > Q (3.7) 

then similarly to (13.61) we get 

+y^ +u^ + v^ n) = 4, x( n ^+y( n> = 2 and u^+v^ = 2, for any n > 0. (3.8) 
To complete the proof we need to the following 
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Lemma 4. If conditions Ii3.5\) and (3.7) are satisfied then 

y(k) _ y (k) _ Q j Qr £ > Q. 


Proof. From (|2.2I1 we get 


y(k+l) _ l_ x (k) v (k) _|_ v (k+l) 


(3.9) 


(3.10) 


y( k+1 ) = ■ 

Now using (13.81) and (13.71) from (|3.9I) we get 

y(k+ 1 ) = 1(2 — y) k ))y) k ) -j- ^(fc+ 1 ) —: _|_ y)k+ 1 ) 

u( fc+1 ) = (i(2 — t/ fc )) + iy( fc )) y^ = \y^■ 

If y(°) = 0 then from the first equation of (I3.10P we get = 0. Consequently 

the second equation gives v ^ = 0. Then using the first equation we get y ( 2) = 0 and 
so on, we get y^ = 0 for any k > 0 . 

If y(°) = 0 then the second equation gives r/ 1 ) = 0. Assume i/°) = v 0 then from 
the first equation we get y W = v + r/ 1 ) = v. Then v Consequently, y® = ^ v. 
Now condition yl 2 )?/ 2 ) = 0 gives v = 0. This completes the proof. □ 


Now by Lemma [1] and property (13.81) we get 

x = 2 and u^ = 2 , for any n > 0 . 

This completes the proof Lemma [3j □ 

Lemma 5. If s = (x, y, u,v) G P is an initial point with x + y + u + v> 4, for which 

(a) if there exists k> 0 such that {x^ k ) +y^){u^ +v^) < 4 then \3.3\) is satisfied. 

(b) i/max{f, if, > 1 then 

lim W n (s) = oo, i.e. at least one coordinate of W n (s) goes to oo. 

n—> oo 

Proof, (a) This simply follows from the equality 

x( fe+1 ) + y (fc+1) + yf k+V) + u (fc+1) = (x (fc) + y (fc) )(w (fc) + u (fc) ). 

Indeed, from this equality it follows that IF(s^ fc+1 ^) G Q 4 . Since Q 4 is invariant the 
part (a) follows from Lemma [2j 

(b) Let us prove it for the case max{ “. qf, = x > 1- For other cases the proof 

is similar. From (12.21) for any s = (x, y , u,v) G P we get 

x [k+1) > ^x (k) u {k \ w (fc+1) > k> 0. (3.11) 

Iterating these inequalities we obtain 

x (fc+1) > ^x (fe) u (fe) > 2 -( 1 + 2 + 2 2 +-+ 2 '=)( xu ) 2fc = 2 , k > 0. (3.12) 

Similarly 

u ( fc + i ) > 2 (^) 2 , k > 0 . 

This completes the proof. □ 
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Denote 

Po = {s = ( X , y,u,v) <£ P : {x + y){u + v) < 4}, 

F = {s = (x,y,u,v) € P : x + y + u + v > 4, max{“, > 1 }. 

Summarizing above-mentioned results we get the following 

Theorem 1. If s = (x,y,u,v) € M 4 is such that 

(i) one of the following conditions is satisfied 

1) s € P 0 ; 

2) s €Q 4 and the condition of part i) of Lemma{3\ is hold; 

3) s £ M, W(s) € P 0 ; 

4) sgU 0 , W 2 (s)gP 0 ; 

5) seJV b W 2 {s) £ P 0 
then 

lirn W n (s) = (0,0, 0,0). 

n—too 

(ii) one of the following conditions is satisfied 

a) s G F; 

b) s G AT, W{s ) € F; 

c) s € No, W 2 (s ) € F; 

d) s G A/i, W 2 (s) € F 
then 

lim lF n '(s) = +oo. 

n—>■ oo 

Proof, (i) The case 1) follows from Lemma [2] and Lemma [5j The case 2) is result 
of Lemma [3j By Lemma CO we have W(M) C F, W 2 (Mo) C F and W 2 (N\) C F. 
Consequently, parts 3)-5) follow from 1 ). 

Part (ii) is a result of Lemma [5] and Lemma [1] 

□ 

Remark 1. The sum of sets for initial points considered in Theorem. [7] is not equal 
to M 4 . But for each fixed point this theorem already gives a large set for the initial 
point, trajectory of which converges to the fixed point. Each point s £ M n can be 
considered as a state of the system, which is a generalized measure (or charge) on the 
set {XX, XX n ,XY, X h Y}. We considered such measure, because for certain purposes, 
it is useful to have a ’’measure” whose values are not restricted to the non-negative reals 
or infinity. Moreover, the dynamical systems considered in this section are interesting 
because they are higher dimensional and such dynamical systems are important, but 
there are relatively few dynamical phenomena that are currently understood JZF. US- 
In the next section we reduce our operators to the invariant set of vectors with non¬ 
negative coordinates (the usual measures which take non-negative values). Then to get 
a stochastic system of probability measures we use a normalization of the non-negative 


measures. 
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4. A NORMALIZED GONOSOMAL OPERTOR 


We note that the gonosomal operator (12.41) does not map the simplex 


S n+U 1 = <s = (x 1 ,...,x n ,y 1 ,...,y v )€ M n+I " : X* > 0, yj > 0, ^ Xj + ^ yj = 1 
[ i =i j =i 

to itself, since 


• Y,y'> = Y 2%' 


(4.1) 


i =1 j=l \ z—1 / \ j = l 

is not equal to 1 in general. 

We denote 

O = {s G S n+U ~ l : (xi,... ,x n ) = (0,.. • ,0)or (y u ... ,y v ) = (0,... ,0)} . 

_ gn-\-v —1 \ ^ 

It is easy to see that W(0) = {(0,... ,0)}. So the points from O do not give any 
contribution to the dynamical system generated by W. 

Therefore we introduce the normalized gonasomal operator as the following. Con¬ 
sider the coefficients of the operator (12.41) with the following properties 

7&>0, 7‘Cl >0, 

+ ELitSj = !. fora]1 

An normalized evolution operator V, with coefficients (|4.2I) is defined as 

Tr-yTl,!; (/) 


(4.2) 


V : 


X 3 = 


y\ = 


££=!*<) 

E n,^ (m) 

»,Ui 7jfc,/gjyfc 

(E”„ **) (EJ.1 w 


j = l,...,n 


, Z = 1,...,p. 


(4.3) 


Proposition 3. The operator V defined by 'with coefficients ^o,ps S n ' v to 

itself if and only if the following condition 


( 7 $,..., 7 £,, 7 $,..., 7 &>)eS”' for all i,k. (4.4) 


is satisfied. 


Proof. Necessity. Suppose for any s € S n ' u we have s' = V(s) € S n,u then we shall 
show that (|4.4p is satisfied. Assume that (|4.4|) is not true, then there is iq £ { 1 ,..., n} 
and ko € { 1 ,... , v} such that 


(lill 0,1 > • • 

• ’ Ti 0 k 0 ,nj 

= ( 0 ,- - 

.., 0 ), 

(4.5) 

l i m ) 

l Uofco,!’ 

(m) \ 

■’7iok o,u) 

= ( 0 ,-. 

■ • , 0 ). 

(4.6) 
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Consider the case (14.5|) (the case (14.61) is similar). Take now some s £ S n,u such that 
Xi 0 / 0, Xi = 0 for i 7 ^ iq and y ko / 0, yk = 0 for k / Uq. Then for this s we have 



= r i 0 k 0 jXi 0 y k0 

(EL, *) (E u ft) (E”=1 x.) (E u ft) 


= 0 , for all j = 1,, n, 


i.e., s' £ O. This is contradiction to the assumption that s' £ S n,v . 

Sufficiency. Assume the conditions (|4.2I) and (14.41) are satisfied, we want to show 
that if s £ S n ' v then s' = U(s) £ S n,L/ . By the construction of the operator (14.31) it is 
easy to see that s' £ 5' n + i/ + 1 so it remains to show that s' O. Assume that s' £ O, i.e, 
(x\ ,..., x' n ) = (0,... , 0) (the case (y[,... , y' u ) = (0,..., 0) is similar). Then by (14.31) 
we should have 

72 , 1 / 

for each 3 = !) ■ ■ •, n. (4.7) 

i,k=1 

Since s £ S n,v there is io £ {l,...,n} and k$ £ such that Xi 0 > 0 and 

( f) 

y ko > 0. From our conditions it follows that 7 i k jXiy k > 0, for all i,j,k. Hence from 
(HU) we get 

( f) 

lioko^ioVko = °> for each 7 = 1,---,^, 

consequently, 

( f) 

Tiokod = 0 for each j = 

This is contradiction to the condition (14.41) . □ 


A fixed point s = (xi,... ,x n ,y\,... ,y v ) of the gonosomal operator (12.41) is called 
non-negative and normalizeable if all coordinates of this point are non-negative and 

£"= i x i + ELi yk > o- 

Proposition 4. There is one-to-one correspondence between non-negative and normal¬ 
izeable fixed points of | and all fixed points of \4-3\) . 

Proof. Let s = (xi,..., x n , y \,..., y u ) be a non-negative and normalizeable fixed point 
of (12.41) . Denote Z = Ya=i + £fc=i 2/fe> and consider the point 

s = {xi/Z ,..., x n /Z , yi/Z ,..., y v /Z). 

By (14.11) for the fixed point we have 

X>Y (4-8) 

k =1 ) 


z = 


u=i j 


Using formula (14.81) one can see that s is a fixed point of (14.31) . 
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Now let s = (x\,... ,x n ,yi,... ,y u ) be a fixed point of (14.31) . i.e. it satisfies the 
following system 


Xj = 


Vi = 


(EM<) (EM* 

E n.v (m) ~ ~ 

i,k=\Tik,l x iyk 


Denote 


(EM<) (EJ-i n 

/ n \ / v 

EM IE 


3 = 1, 


Z = l, 


,n 


(4.9) 


Z = 


Vk 


(4.10) 


U=i 


vfc=l 


Dividing both side of (14.91) to Z it is easy to see that the following point is a fixed point 
of {21} : s = ( 'xi/Z ,.. ,,x n /Z,yi/Z ,.. .,y v /Z). □ 

In this section we consider the normalized version of the evolution operator (12.21) . 
i.e., 

2 xu + yu 


V : 


x = 


y = 


u = 


V = 


4(x + y)(it + v) 

6am + 3 yu + 4 yv 
12(x + y)(w + v) 

6 xu + 6am + 3 yu + 4 yv 
12(x + y)(rt + v ) 
2>yu + 4yv 


(4.11) 


12(x + y)(w + v)' 

It is easy to see that the operator (14.111) satisfies the conditions of Proposition [3l hence 
V : S 2 ’ 2 -A- S 2 ’ 2 . 

The following lemmas give some useful estimates. 

Lemma 6. 

14-11 1 ) then 


to 

II 

IT 

e 

) 6 5 2,2 and 


= (s', 

y', 

u', ' 


U 

< x' 

< 

U 

< 

1 


4 (u+v) 

2(u-\-v) 

2’ 


V 

< y' 

< 

u-\-2v 

< 

1 


3 {u+v) 

4(u-\-v) 

2 ’ 

i < 

2x+y 

4 (x+y) 

< v! 

< 

3x+2y 

6(x+y) 

< 

1 

2 ’ 


y 

< i/ 

< 

y 

< 

1 


4 (x+y) 

3 (x+y) 

3’ 

5 + 

u 

< x' + y' 

< 

1 



6 (u+v) 

2’ 




^ < v! + v' < ^ + 

/ 


2 — 
v' < 




< I. 


— 2 1 — 3’ 

yr < u', 

x' < u'. 
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Proof. Straightforward. 

Lemma 7. Let s = (x,y,u,v) £ <S 2,2 and = {x^ n \y^ n \u^ n \ 
operator then 

1. 

A < X (2) ly(2) < I. 

12 ~ y - 2’ 

2. There exists a £ (0,1) such that 

p( n+1 ) < ay( n \ n > 2. 


Proof. 1. Using Lemma [6] we have 

+ 


— H—-———- < x^ + y^ < —. 
3 6 {u' + v')~ y ~ 2 


Consequently, since u! > i/ we get 

1 

- + 


> 1 1 - 5 


3 6(ii' + U) _ 3 12 12' 

2. Using u = l — x — y — v we can rewrite the operator (I4.11|) 


V : 


l (2x+y)(l-x-y-v) 

4(x+y)(l-x-y) 

I _ 3y(l-x-y)+(6x+y)v 
■' 12(x+y)(l-x-y) 

/ 3(1— x—y)+v 

V = V 12(l + y)(l-x-y) - 


For any n > 1 we have 


V 


71+1 


where 
<p(x,y) = 


x 


' (n+l) = (2xW+y("))(l-^M-yW- t ;W) 

4(a?( Tl )+2/( n ))(l— x( n ^ —y ( n )) 

y(n+l) _ v (n+l) _|_ ipfx^ n \ y n )v^ 
^In+i) _ ijj( x ( n ) ^ y( n ))y( n ) ^ 

3(1 — a: — y) + v 




2(x+ y)(l - x - y) \2{x + y){l - x - y) 

Using above mentioned inequalities we get 

1 ri 

+£,?/) <U ip{x,v) = 


3(1 —u — v) 12(zi + u)(l — u — v) 
Consequently we get from (14.131) the following 

{ y(n+ 1) < v (n+l) _|_ v (n) 
v (n+l) < 13 y (n)_ 


□ 

+)) = V n (s ) /or t/ie 


(4.12) 

(4.13) 

3n + Av 

12(w + u)(l — u — v) 

2 1 _ 13 

“ 3 ~~ 8 “ 24' 

(4.14) 


This completes the proof. 


□ 
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Theorem 2. The operator \f.ll\ ) has a unique fixed point p = (1/2,0,1/2, 0) and there 
is an open neighborhood U(p) C S 2,2 of p such that for any initial point s £ U{jp) we 
have 


lim U n (s) = p. 

n—>■ oo 


Proof. The existence and uniqueness of p follow from Propositions Q] and [4J To check 
the Jacobi matrix of the operator (14.111) at the fixed point p, one has to replace v by 
v = 1 — x — y — it in the operator then construct a 3 x 3 Jacobi matrix. It is easy to see 
that this matrix at the fixed point has eigenvalues —0.5, 0,1, i.e the point is attractive, 
so HJ Theorem 6.3] completes the proof. □ 


Using Lemma [6] and Lemma 0 one can see that the trajectory of any initial point 
after few iterations comes close to the fixed point p. Moreover, using Maple one can 
see that the limit point of the trajectory is always p. Thus the following should be true 


Conjecture. For any initial point s £ S 2,2 we have lim n _ > , 00 V n (s) = p. 

Remark 2. The results have the following biological interpretations: Let s = (x, y, u, v ) £ 
S 2 ' 2 be an initial state (the probability distribution on the set {XX,XX h ;XY,X h Y} 
of genotypes). Theorem^ says that, as a rule, the population tends to the equilibrium 
state p = (1/2,0,1/2,0) with the passage of time, i.e. the future of the population is 
stable: genotypes XX and XY are survived always, but the genotypes XX h and X h Y 
(therefore hemophilia) will disappear in the future. It follows that hemophilia is main¬ 
tained in a population only if it occurs mutations on the genes coding for the coagulation 
factors. 
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